Complete mitochondrial genome sequencing reveals novel haplotypes in a Polynesian population

The high risk of metabolic disease traits in Polynesians may be partly explained by elevated prevalence of genetic variants involved in energy metabolism. The genetics of Polynesian populations has been shaped by island hoping migration events which have possibly favoured thrifty genes. The aim of this study was to sequence the mitochondrial genome in a group of Maoris in an effort to characterise genome variation in this Polynesian population for use in future disease association studies. We sequenced the complete mitochondrial genomes of 20 non-admixed Maori subjects using Affymetrix technology. DNA diversity analyses showed the Maori group exhibited reduced mitochondrial genome diversity compared to other worldwide populations, which is consistent with historical bottleneck and founder effects. Global phylogenetic analysis positioned these Maori subjects specifically within mitochondrial haplogroup - B4a1a1. Interestingly, we identified several novel variants that collectively form new and unique Maori motifs – B4a1a1c, B4a1a1a3 and B4a1a1a5. Compared to ancestral populations we observed an increased frequency of non-synonymous coding variants of several mitochondrial genes in the Maori group, which may be a result of positive selection and/or genetic drift effects. In conclusion, this study reports the first complete mitochondrial genome sequence data for a Maori population. Overall, these new data reveal novel mitochondrial genome signatures in this Polynesian population and enhance the phylogenetic picture of maternal ancestry in Oceania. The increased frequency of several mitochondrial coding variants makes them good candidates for future studies aimed at assessment of metabolic disease risk in Polynesian populations.

Matters Arising: HLA-DR-DQ haplotypes and genotypes in Finnish patients with rheumatoid arthritis

We read with great interest the paper by Laivoranta-Nyman et al.1 They studied Finnish patients with rheumatoid arthritis (RA) and controls, and suggested that there exist susceptibility, neutral and protective HLA-DR-DQ haplotypes that do not match the predictions of current hypotheses for the mechanism of association of the HLA class II region and RA...

Dissection of class III major histocompatibility complex haplotypes associated with rheumatoid arthritis

Objective. We have previously identified a single-nucleotide polymorphism (SNP) haplotype involving the lymphotoxin α (LTA) and tumor necrosis factor (TNF) loci (termed haplotype LTA-TNF2) on chromosome 6 that shows differential association with rheumatoid arthritis (RA) on HLA-DRB1*0404 and *0401 haplotypes, suggesting the presence of additional non-HLA-DRB1 RA susceptibility genes on these haplotypes. To refine this association, we performed a case-control association study using both SNPs and microsatellite markers in haplotypes matched either for HLA-DRB1*0404 or for HLA-DRB1*0401. Methods. Fourteen SNPs lying between HLA-DRB1 and LTA were genotyped in 87 DRB1*04-positive families. High-density microsatellite typing was performed using 24 markers spanning 2,500 kb centered around the TNF gene in 305 DRB1*0401 or *0404 cases and 400 DRB1*0401 or *0404 controls. Single-marker, 2-marker, and 3-marker minihaplotypes were constructed and their frequencies compared between the DRB1*0401 and DRB1*0404 matched case and control haplotypes. Results. Marked preservation of major histocompatibility complex haplotypes was seen, with chromosomes carrying LTA-TNF2 and either DRB1*0401 or DRB1*0404 both carrying an identical SNP haplotype across the 1-Mb region between TNF and HLA-DRB1. Using microsatellite markers, we observed two 3-marker minihaplotypes that were significantly overrepresented in the DRB1*0404 case haplotypes (P = 0.00024 and P = 0.00097). Conclusion. The presence of a single extended SNP haplotype between LTA-TNF2 and both DRB1*0401 and DRB1*0404 is evidence against this region harboring the genetic effects in linkage disequillbrium with LTA-TNF2. Two RA-associated haplotypes on the background of DRB1*0404 were identified in a 126-kb region surrounding and centromeric to the TNF locus.

Letter Reply: Distribution of TNFA haplotypes in healthy Caucasians: Comment on the articles by Newton et al and Zeggini et al

We thank Ploski and colleagues for their interest in our study. The explanation for the difference in our findings is a typographic error in Table 2 of our article, whereby the alleles for marker TNF ⫺1031 were labeled incorrectly...

Dissection of the genetics of Parkinson's disease identifies an additional association 5' of SNCA and multiple associated haplotypes at 17q21

We performed a genome-wide association study (GWAS) in 1705 Parkinson's disease (PD) UK patients and 5175 UK controls, the largest sample size so far for a PD GWAS. Replication was attempted in an additional cohort of 1039 French PD cases and 1984 controls for the 27 regions showing the strongest evidence of association (P < 10 4). We replicated published associations in the 4q22/SNCA and 17q21/MAPT chromosome regions (P < 10 10) and found evidence for an additional independent association in 4q22/SNCA.A detailed analysis of the haplotype structure at 17q21 showed that there are three separate risk groups within this region. We found weak but consistent evidence of association for common variants located in three previously published associated regions (4p15/BST1, 4p16/GAK and 1q32/PARK16). We found no support for the previously reported SNP association in 12q12/LRRK2. We also found an association of the two SNPs in 4q22/SNCA with the age of onset of the disease. © The Author 2010. Published by Oxford University Press.

Serum bilirubin concentration is modified by UGT1A1 Haplotypes and influences risk of type-2 diabetes in the Norfolk Island Genetic Isolate

Background Located in the Pacific Ocean between Australia and New Zealand, the unique population isolate of Norfolk Island has been shown to exhibit increased prevalence of metabolic disorders (type-2 diabetes, cardiovascular disease) compared to mainland Australia. We investigated this well-established genetic isolate, utilising its unique genomic structure to increase the ability to detect related genetic markers. A pedigree-based genome-wide association study of 16 routinely collected blood-based clinical traits in 382 Norfolk Island individuals was performed. Results A striking association peak was located at chromosome 2q37.1 for both total bilirubin and direct bilirubin, with 29 SNPs reaching statistical significance (P < 1.84 × 10−7). Strong linkage disequilibrium was observed across a 200 kb region spanning the UDP-glucuronosyltransferase family, including UGT1A1, an enzyme known to metabolise bilirubin. Given the epidemiological literature suggesting negative association between CVD-risk and serum bilirubin we further explored potential associations using stepwise multivariate regression, revealing significant association between direct bilirubin concentration and type-2 diabetes risk. In the Norfolk Island cohort increased direct bilirubin was associated with a 28 % reduction in type-2 diabetes risk (OR: 0.72, 95 % CI: 0.57-0.91, P = 0.005). When adjusted for genotypic effects the overall model was validated, with the adjusted model predicting a 30 % reduction in type-2 diabetes risk with increasing direct bilirubin concentrations (OR: 0.70, 95 % CI: 0.53-0.89, P = 0.0001). Conclusions In summary, a pedigree-based GWAS of blood-based clinical traits in the Norfolk Island population has identified variants within the UDPGT family directly associated with serum bilirubin levels, which is in turn implicated with reduced risk of developing type-2 diabetes within this population.

Genetic resources in wild and cultured stocks of the Asian mud crab, Scylla paramamosain

The mud crab (Scylla spp.) aquaculture industry has expanded rapidly in recent years in many countries in the Indo - West Pacific (IWP) region as an alternative to marine shrimp culture because of significant disease outbreaks and associated failures of many shrimp culture industries in the region. Currently, practices used to produce and manage breeding crabs in hatcheries may compromise levels of genetic diversity, ultimately compromising growth rates, disease resistance and stock productivity. Therefore, to avoid “genetic pollution” and its harmful effects and to promote further development of mud crab aquaculture and fisheries in a sustainable way, a greater understanding of the genetic attributes of wild and cultured mud crab stocks is required. Application of these results can provide benefits for managing wild and cultured Asian mud crab populations for multiple purposes including for commercial production, recreation and conservation and to increase profitability and sustainability of newly emerging crab culture industries. Phylogeographic patterns and the genetic structure of Asian mud crab populations across the IWP were assessed to determine if they were concordant with those of other widespread taxa possessing pelagic larvae of relatively long duration. A 597 bp fragment of the mitochondrial DNA COI gene was amplified and screened for variation in a total of 297 individuals of S. paramamosain from six sampling sites across the species’ natural geographical distribution in the IWP and 36 unique haplotypes were identified. Haplotype diversities per site ranged from 0.516 to 0.879. Nucleotide diversity estimates among haplotypes were 0.11% – 0.48%. Maximum divergence observed among S. paramamosain samples was 1.533% and samples formed essentially a single monophyletic group as no obvious clades were related to geographical location of sites. A weak positive relationship was observed however, between genetic distance and geographical distance among sites. Microsatellite markers were then used to assess contemporary gene flow and population structure in Asian mud crab populations sampled across their natural distribution in the IWP. Eight microsatellite loci were screened in sampled S. paramamosain populations and all showed high allelic diversity at all loci in sampled populations. In total, 344 individuals were analysed, and 304 microsatellite alleles were found across the 8 loci. The mean number of alleles per locus at each site ranged from 20.75 to 28.25. Mean allelic richness per site varied from 17.2 to 18.9. All sites showed high levels of heterozygosity as average expected heterozygosities for all loci ranged from 0.917 – 0.953 while mean observed heterozygosity ranged from 0.916 – 0.959. Allele diversities were similar at all sites and across all loci. The results did not show any evidence for major differences in allele frequencies among sites and patterns of allele frequencies were very similar in all populations across all loci. Estimates of population differentiation (FST) were relatively low and most probably largely reflect intra – individual variation for very highly variable loci. Results from nDNA analysis showed evidence for only very limited population genetic structure among sampled S. paramamosain, and a positive and significant association for genetic and geographical distance among sample sites. Microsatellite markers were then employed to determine if adequate levels of genetic diversity has been captured in crab hatcheries for the breeding cycle. The results showed that all microsatellite loci were polymorphic in hatchery samples. Culture populations were in general, highly genetically depauperate, compared with comparable wild populations, with only 3 to 8 alleles recorded for the same loci set per population. In contrast, very high numbers of alleles per locus were found in reference wild S. paramamosain populations, which ranged from 18 to 46 alleles per locus per population. In general, this translates into a 3 to 10 fold decline in mean allelic richness per locus in all culture stocks compared with wild reference counterparts. Furthermore, most loci in all cultured S. paramamosain samples showed departures from HWE equilibrium. Allele frequencies were very different in culture samples from that present in comparable wild reference samples and this in particular, was reflected in a large decline in allele diversity per locus. The pattern observed was best explained by significant impacts of breeding practices employed in hatcheries rather than natural differentiation among wild populations used as the source of brood stock. Recognition of current problems and management strategies for the species both for the medium and long-term development of the new culture industry are discussed. The priority research to be undertaken over the medium term for S. paramamosain should be to close the life cycle fully to allow individuals to be bred on demand and their offspring equalised to control broodstock reproductive contributions. Establishing a broodstock register and pedigree mating system will be required before any selection program is implemented. This will ensure that sufficient genetic variation will be available to allow genetic gains to be sustainably achieved in a future stock improvement program. A fundamental starting point to improve hatchery practices will be to encourage farmers and hatchery managers to spawn more females in their hatcheries as it will increase background genetic diversity in culture stocks. Combining crablet cohorts from multiple hatcheries into a single cohort for supply to farmers or rotation of breeding females regularly in hatcheries will help to address immediate genetic diversity problems in culture stocks. Application of these results can provide benefits for managing wild and cultured Asian mud crab populations more efficiently. Over the long-term, application of data on genetic diversity in wild and cultured stocks of Asian mud crab will contribute to development of sustainable and productive culture industries in Vietnam and other countries in the IWP and can contribute towards conservation of wild genetic resources.

Dysbindin (DTNBP1) : a role in psychotic depression?

Previous studies yielded evidence for dysbindin (DTNBP1) to impact the pathogenesis of schizophrenia on the one hand and affective disorders such as bipolar or major depressive disorder (MDD) on the other. Thus, in the present study we investigated whether DTNBP1 variation was associated with psychotic depression as a severe clinical manifestation of MDD possibly constituting an overlapping phenotype between affective disorders and schizophrenia. A sample of 243 Caucasian inpatients with MDD (SCID-I) was genotyped for 12 SNPs spanning 92% of the DTNBP1 gene region. Differences in DTNBP1 genotype distributions across diagnostic subgroups of psychotic (N = 131) vs. non-psychotic depression were estimated by Pearson Chi2 test and logistic regression analyses adjusted for age, gender, Beck Depression Inventory (BDI) and the Global Assessment of Functioning Scale (GAF). Overall, patients with psychotic depression presented with higher BDI and lower GAF scores expressing a higher severity of the illness as compared to depressed patients without psychotic features. Four DTNBP1 SNPs, particularly rs1997679 and rs9370822, and the corresponding haplotypes, respectively, were found to be significantly associated with the risk of psychotic depression in an allele-dose fashion. In summary, the present results provide preliminary support for dysbindin (DTNBP1) gene variation, particularly SNPs rs1997679 and rs9370822, to be associated with the clinical phenotype of psychotic depression suggesting a possible neurobiological mechanism for an intermediate trait on the continuum between affective disorders and schizophrenia.

Evolutionary biology of Gondwanan non-biting midges (Diptera: Chironomidae)

The potential restriction to effective dispersal and gene flow caused by habitat fragmentation can apply to multiple levels of evolutionary scale; from the fragmentation of ancient supercontinents driving diversification and speciation on disjunct landmasses, to the isolation of proximate populations as a result of their inability to cross intervening unsuitable habitat. Investigating the role of habitat fragmentation in driving diversity within and among taxa can thus include inferences of phylogenetic relationships among taxa, assessments of intraspecific phylogeographic structure and analyses of gene flow among neighbouring populations. The proposed Gondwanan clade within the chironomid (non-biting midge) subfamily Orthocladiinae (Diptera: Chironomidae) represents a model system for investigating the role that population fragmentation and isolation has played at different evolutionary scales. A pilot study by Krosch et al (2009) indentified several highly divergent lineages restricted to ancient rainforest refugia and limited gene flow among proximate sites within a refuge for one member of this clade, Echinocladius martini Cranston. This study provided a framework for investigating the evolutionary history of this taxon and its relatives more thoroughly. Populations of E. martini were sampled in the Paluma bioregion of northeast Queensland to investigate patterns of fine-scale within- and among-stream dispersal and gene flow within a refuge more rigorously. Data was incorporated from Krosch et al (2009) and additional sites were sampled up- and downstream of the original sites. Analyses of genetic structure revealed strong natal site fidelity and high genetic structure among geographically proximate streams. Little evidence was found for regular headwater exchange among upstream sites, but there was distinct evidence for rare adult flight among sites on separate stream reaches. Overall, however, the distribution of shared haplotypes implied that both larval and adult dispersal was largely limited to the natal stream channel. Patterns of regional phylogeographic structure were examined in two related austral orthoclad taxa – Naonella forsythi Boothroyd from New Zealand and Ferringtonia patagonica Sæther and Andersen from southern South America – to provide a comparison with patterns revealed in their close relative E. martini. Both taxa inhabit tectonically active areas of the southern hemisphere that have also experienced several glaciation events throughout the Plio-Pleistocene that are thought to have affected population structure dramatically in many taxa. Four highly divergent lineages estimated to have diverged since the late Miocene were revealed in each taxon, mirroring patterns in E. martini; however, there was no evidence for local geographical endemism, implying substantial range expansion post-diversification. The differences in pattern evident among the three related taxa were suggested to have been influenced by variation in the responses of closed forest habitat to climatic fluctuations during interglacial periods across the three landmasses. Phylogeographic structure in E. martini was resolved at a continental scale by expanding upon the sampling design of Krosch et al (2009) to encompass populations in southeast Queensland, New South Wales and Victoria. Patterns of phylogeographic structure were consistent with expectations and several previously unrecognised lineages were revealed from central- and southern Australia that were geographically endemic to closed forest refugia. Estimated divergence times were congruent with the timing of Plio-Pleistocene rainforest contractions across the east coast of Australia. This suggested that dispersal and gene flow of E. martini among isolated refugia was highly restricted and that this taxon was susceptible to the impacts of habitat change. Broader phylogenetic relationships among taxa considered to be members of this Gondwanan orthoclad group were resolved in order to test expected patterns of evolutionary affinities across the austral continents. The inferred phylogeny and estimated divergence times did not accord with expected patterns based on the geological sequence of break-up of the Gondwanan supercontinent and implied instead several transoceanic dispersal events post-vicariance. Difficulties in appropriate taxonomic sampling and accurate calibration of molecular phylogenies notwithstanding, the sampling regime implemented in the current study has been the most intensive yet performed for austral members of the Orthocladiinae and unsurprisingly has revealed both novel taxa and phylogenetic relationships within and among described genera. Several novel associations between life stages are made here for both described and previously unknown taxa. Investigating evolutionary relationships within and among members of this clade of proposed Gondwanan orthoclad taxa has demonstrated that a complex interaction between historical population fragmentation and dispersal at several levels of evolutionary scale has been important in driving diversification in this group. While interruptions to migration, colonisation and gene flow driven by population fragmentation have clearly contributed to the development and maintenance of much of the diversity present in this group, long-distance dispersal has also played a role in influencing diversification of continental biotas and facilitating gene flow among disjunct populations.

The factor V HR2 haplotype: Prevalence and association of the A4070G and A6755G polymorphisms

Recently, a polymorphism was identified in exon 25 of the factor V gene that is possibly a functional candidate for the HR2 haplotype. This haplotype is characterized by a single base substitution named R2 (A4070G) in the B domain of the protein. A mutation (A6755G; 2194Asp→Gly) located near the C terminus has been hypothesized to influence protein folding and glycosylation, and might be responsible for the shift in factor V isoform (FV1 / FV2) ratio. This study investigated the prevalence of these two factor V HR2 haplotype polymorphisms in a cohort of normal blood donors, patients with osteoarthritis and women with complications during pregnancy, and in families of factor V Leiden individuals. A high allele frequency for the two polymorphisms was found in the blood donor group (6.2% R2, 5.6% A6755G). No significant difference in allele frequency was observed in the clinical groups (obstetric complications and osteoarthritis, 4.1-4.9% for the two polymorphisms) when compared with that of healthy blood donors. We confirm that the factor V A6755G polymorphism shows strong linkage to the R2 allele, although it is not exclusively inherited with the exon 13 A4070G variant and can occur independently. © 2001 Lippincott Williams & Wilkins.

Common variation in the fibroblast growth factor receptor 2 gene is not associated with endometriosis risk

BACKGROUND Endometriosis is a polygenic disease with a complex and multifactorial aetiology that affects 8-10% of women of reproductive age. Epidemiological data support a link between endometriosis and cancers of the reproductive tract. Fibroblast growth factor receptor 2 (FGFR2) has recently been implicated in both endometrial and breast cancer. Our previous studies on endometriosis identified significant linkage to a novel susceptibility locus on chromosome 10q26 and the FGFR2 gene maps within this linkage region. We therefore hypothesized that variation in FGFR2 may contribute to the risk of endometriosis. METHODS We genotyped 13 single nucleotide polymorphisms (SNPs) densely covering a 27 kb region within intron 2 of FGFR2 including two SNPs (rs2981582 and rs1219648) significantly associated with breast cancer and a total 40 tagSNPs across 150 kb of the FGFR2 gene. SNPs were genotyped in 958 endometriosis cases and 959 unrelated controls. RESULTS We found no evidence for association between endometriosis and FGFR2 intron 2 SNPs or SNP haplotypes and no evidence for association between endometriosis and variation across the FGFR2 gene. CONCLUSIONS Common variation in the breast-cancer implicated intron 2 and other highly plausible causative candidate regions of FGFR2 do not appear to be a major contributor to endometriosis susceptibility in our large Australian sample.

Analysis of the CDKN2A, CDKN2B and CDK4 genes in 48 Australian melanoma kindreds

Germline mutations within the cyclin-dependent kinase inhibitor 2A (CDKN2A) gene and one of its targets, the cyclin dependent kinase 4 (CDK4) gene, have been identified in a proportion of melanoma kindreds. In the case of CDK4, only one specific mutation, resulting in the substitution of a cysteine for an arginine at codon 24 (R24C), has been found to be associated with melanoma. We have previously reported the identification of germline CDKN2A mutations in 7/18 Australian melanoma kindreds and the absence of the R24C CDK4 mutation in 21 families lacking evidence of a CDKN2A mutation. The current study represents an expansion of these efforts and includes a total of 48 melanoma families from Australia. All of these families have now been screened for mutations within CDKN2A and CDK4, as well as for mutations within the CDKN2A homolog and 9p21 neighbor, the CDKN2B gene, and the alternative exon 1 (E1beta) of CDKN2A. Families lacking CDKN2A mutations, but positive for a polymorphism(s) within this gene, were further evaluated to determine if their disease was associated with transcriptional silencing of one CDKN2A allele. Overall, CDKN2A mutations were detected in 3/30 (10%) of the new kindreds. Two of these mutations have been observed previously: a 24 bp duplication at the 5' end of the gene and a G to C transversion in exon 2 resulting in an M531 substitution. A novel G to A transition in exon 2, resulting in a D108N substitution was also detected. Combined with our previous findings, we have now detected germline CDKN2A mutations in 10/48 (21%) of our melanoma kindreds. In none of the 'CDKN2A-negative' families was melanoma found to segregate with either an untranscribed CDKN2A allele, an R24C CDK4 mutation, a CDKN2B mutation, or an E1beta mutation. The last three observations suggest that these other cell cycle control genes (or alternative gene products) are either not involved at all, or to any great extent, in melanoma predisposition.

Haplotype analysis of two recurrent CDKN2A mutations in 10 melanoma families : evidence for common founders and independent mutations

Germ-line mutations in CDKN2A have been shown to predispose to cutaneous malignant melanoma. We have identified 2 new melanoma kindreds which carry a duplication of a 24bp repeat present in the 5' region of CDKN2A previously identified in melanoma families from Australia and the United States. This mutation has now been reported in 5 melanoma families from 3 continents: Europe, North America, and Australasia. The M53I mutation in exon 2 of CDKN2A has also been documented in 5 melanoma families from Australia and North America. The aim of this study was to determine whether the occurrence of the mutations in these families from geographically diverse populations represented mutation hotspots within CDKN2A or were due to common ancestors. Haplotypes of 11 microsatellite markers flanking CDKN2A were constructed in 5 families carrying the M53I mutation and 5 families carrying the 24bp duplication. There were some differences in the segregating haplotypes due primarily to recombinations and mutations within the short tandem-repeat markers; however, the data provide evidence to indicate that there were at least 3 independent 24bp duplication events and possibly only 1 original M53I mutation. This is the first study to date which indicates common founders in melanoma families from different continents.

The ghrelin signalling system is involved in the consumption of sweets

The gastric-derived orexigenic peptide ghrelin affects brain circuits involved in energy balance as well as in reward. Indeed, ghrelin activates an important reward circuit involved in natural- as well as drug-induced reward, the cholinergic-dopaminergic reward link. It has been hypothesized that there is a common reward mechanism for alcohol and sweet substances in both animals and humans. Alcohol dependent individuals have higher craving for sweets than do healthy controls and the hedonic response to sweet taste may, at least in part, depend on genetic factors. Rat selectively bred for high sucrose intake have higher alcohol consumption than non-sucrose preferring rats and vice versa. In the present study a group of alcohol-consuming individuals selected from a population cohort was investigated for genetic variants of the ghrelin signalling system in relation to both their alcohol and sucrose consumption. Moreover, the effects of GHS-R1A antagonism on voluntary sucrose- intake and operant self-administration, as well as saccharin intake were investigated in preclinical studies using rodents. The effects of peripheral grelin administration on sucrose intake were also examined. Here we found associations with the ghrelin gene haplotypes and increased sucrose consumption, and a trend for the same association was seen in the high alcohol consumers. The preclinical data show that a GHS-R1A antagonist reduces the intake and self-administration of sucrose in rats as well as saccharin intake in mice. Further, ghrelin increases the intake of sucrose in rats. Collectively, our data provide a clear indication that the GHS-R1A antagonists reduces and ghrelin increases the intake of rewarding substances and hence, the central ghrelin signalling system provides a novel target for the development of drug strategies to treat addictive behaviours. © 2011 Landgren et al.

A novel SNP in COMT is associated with alcohol dependence but not opiate or nicotine dependence: a case control study

Background It is well established that COMT is a strong candidate gene for substance use disorder and schizophrenia. Recently we identified two SNPs in COMT (rs4680 and rs165774) that are associated with schizophrenia in an Australian cohort. Individuals with schizophrenia were more than twice as likely to carry the GG genotype compared to the AA genotype for both the rs165774 and rs4680 SNPs. Association of both rs4680 and rs165774 with substance dependence, a common comorbidity of schizophrenia has not been investigated. Methods To determine whether COMT is important in substance dependence, rs165774 and rs4680 were genotyped and haplotyped in patients with nicotine, alcohol and opiate dependence. Results The rs165774 SNP was associated with alcohol dependence. However, it was not associated with nicotine or opiate dependence. Individuals with alcohol dependence were more than twice as likely to carry the GG or AG genotypes compared to the AA genotype, indicating a dominant mode of inheritance. The rs4680 SNP showed a weak association with alcohol dependence at the allele level that did not reach significance at the genotype level but it was not associated with nicotine or opiate dependence. Analysis of rs165774/rs4680 haplotypes also revealed association with alcohol dependence with the G/G haplotype being almost 1.5 times more common in alcohol-dependent cases. Conclusions Our study provides further support for the importance of the COMT in alcohol dependence in addition to schizophrenia. It is possible that the rs165774 SNP, in combination with rs4680, results in a common molecular variant of COMT that contributes to schizophrenia and alcohol dependence susceptibility. This is potentially important for future studies of comorbidity. As our participant numbers are limited our observations should be viewed with caution until they are independently replicated.